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It has been shown that the theory of hadronic tran- 
sitions based on QCD multipole expansion [l], 0, 0, [1] 
can make quite successful predictions for many hadronic 
transition rates in the cc and bb systems 0, So far, 
hadronic transitions of the Xcj states have not been stud- 
ied yet. The observed Xcj decays are mainly hadronic 
decays (decaying into light hadrons), and the hadronic 
widths of the three Xcj states are rather different. From 
the point of view of perturbative QCD, the hadronic de- 
cay rates of XcO and Xc2 are proportional to a1 , while that 
of Xci is proportional to Q ■ Hence Xci has the small- 
est hadronic decay rate, and is thus the most interesting 
one among the three Xcj states for studying hadronic 
transitions. The main hadronic transition process of Xci 
is Xci ~^ ?7c7r7r. In the framework of QCD multipole ex- 
pansion, this is dominated by the El-Ml transition, and 
the rate is significantly smaller than those of El-El tran- 
sitions such as ij:' — > J/il) tttt. In addition, Xci cannot 
be directly produced in e^e'^ collision. It can only be 
produced via the radiative transition ip' — > 7Xci- This 
is why Xci TjcTTTT is not easy to detect at the original 
Beijing Spectrometer BES II. The designed luminosity of 
the upgraded Beijing Electron-Positron Collider BEPC II 
is about 1 — 2 orders of magnitude higher than that of 
the original BEPC I, and the ability of the upgraded Bei- 
jing Spectrometer BES III will be significantly improved 
compared with the original BES II, especially its photon 
detector. So that BES III will be able to detect processes 
like Xci — * 'ycTTTT which BES II can hardly do. In this pa- 
per, we calculate the transition rate of Xci T]cmr in the 
framework of QCD multipole expansion, and our result 
shows that both the upgraded BES III and the Cornell 



CLEO-c have a good chance to study Xci — * 'Hc'K'k. and 
the theoretical approaches to hadronic transitions can be 
further tested. ^ 

According to the formula for El-Ml transitions 2], the 
transition amplitude for Xci Vc''^''^ is 
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where i?" and are color electric and magnetic fields; 
gE and qm are the corresponding coupling constants for 
El and Ml gluon emissions; x is the relative separation 
between c and c in the quarkonium; Sc (ss) is the spin 
of the c (c) quark; Ej is the energy of the initial state 
Xci; Ekl and \KL) are the energy and eigenstate of the 
intermediate state with the principal quantum number 
K and orbital angular momentum L, respectively. As 
described in Ref. we take the quark confining string 
model Q to describe the intermediate states. 

The hadronization factor {TnT\Ef Bj\0) is a second rank 
pseudo tensor at the light hadron scale. In the rest frame 
of Xci, it is a function of the two pion momenta qi and 



(72- Since the mass difference My 



530.3 MeV 



is small, the pions are soft. Taking the PCAC and soft 
pion approach, the hadronization factor is of the form 
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where C is approximately a constant characterizing the 
size of the matrix element. Unfortunately, there is no 
experimental result serving as the input data to deter- 
mine C at present. Therefore we have to take an alter- 
native approach to make predictions. Now we take the 
two gluon approximation used in Refs. 0, which is 
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shown to be reasonable so far as only the transition rates 
are concerned 0,0]. In this approach, we approximately 
express {■KTT\Ef B'^\{)) as 

(w|i?fBf |0) « {nn\gg){gg\EfB'^\Q). (3) 

The matrix element (Trvrlgy) is at the scale of M^^^ — 
Mrj^ ~ 530 MeV, and is independent of the multi- 
pole gluon emission. If we take the same approxima- 
tion to the process ip' —> J/ip tttt, i.e., expressing the 
corresponding hadronization factor as {TnT\E!j^Ej\0) w 
{TTTrlgg) {gg\E^ Ej\0) , the matrix element {'KTi\gg) in this 
process is at the scale of M^i — Mj/^ ~ 590 MeV 
which is not much different from 530 MeV in the case 
of Xci '/cTTTT. Therefore we can neglect the running 
of {■KTT\gg) from 590 MeV to 530 MeV, and regard it 
as the same matrix element in both ij:' —>■ J /ip tttt and 
Xci ~* ?7c7r7r. Then {■m:\gg) can be absorbed into the 



definition of the phenomenological coupling constant gE, 
and we can determine g_E by taking the experimental da- 
tum of the transition rate r(?/'' — > J/ip tttt) as input. 
The remaining factor {gg\EfB'^\Q) in Q is then easy to 
evaluate. It has been shown in Refs. 0, 9 that such an 
approach gives predictions for the transition rates quite 
close to those in the soft pion approach Q. So we take 
this approximation and obtain the following transition 
rate 
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in which i?/, Rp, Rkl are radial wave functions of the 
initial, final, and intermediate states, respectively. These 
radial wave functions can be obtained when a potential 
model for the quarkonium is taken. 

In this paper, we take the improved QCD motivated 
potential model proposed by Chen and Kuang (CK) [TT| . 
and we take Potential I in Ref. Comparing with 

other potential models, this potential model has the ad- 
vantage that it reflects more about QCD with an ex- 
plicit Ajy-j dependence, and leads to more successful phe- 



nomenological results [1, |Tl| . The potential is 
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where k = 0.1491 GeV^ is the string tension related to 
the Regge slope [III, 1e is the Euler constant, and /(r) 
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in which A-^— - = 180 MeV. In this model, the c quark 

MS ^ 

mass is TOc = 1.478 GeV, and we take A-g-g- = 200 MeV 
in the following calculation. 

With this potential model, we can calculate the tran- 
sition ampUtudes fl;fi\^li^{K) defined in Eq. ©. To 

see the convergence of the summation fnfi'jnpiF 

K 

I 



in ([5]), we list the values of the first eight amplitudes in 
TABLE I. We see that taking the first eight terms in the 
summation is enough, so we obtain 

/mo + /mo = -7.21402 GeV-^. (8) 
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TABLE I: The first eight transition amphtudes (in units of GeV obtained from the CK potential model. 





mUK) (GeV-2) 


mUK) (GeV-^) 


fil%{K) + mUK) (GeV-2) 


K=l 


-2.43573 


-5.14597 


-7.58170 


K=2 


-0.17284 


0.50817 


0.33533 


K=3 


-0.02746 


0.05487 


0.02740 


K=4 


-0.00683 


0.01075 


0.00393 


K=5 


-0.00225 


0.00303 


0.00078 


K=6 


-0.00090 


0.00109 


0.00018 


K=7 


-0.00041 


0.00046 


4.5x10"^^ 


K=8 


-0.00021 


0.00022 


9.0x10"'' 


sum 


-2.24676 


-4.56726 


-7.21402 



Next we consider the determination of and um- 
Taking the same approximation to if)' — » J /if) tttt, we can 
obtain T{'ip' ^ J/il) tttt) which is 



r(vy -> J/^ tttt) 



8q; 



\.fllU\M^,~Mj,^y. (9) 



85057r 



The updated Particle Data Group (PDG) best fit exper- 
imental values [13] 

T{tp') = 337 ± 13 keV, 

-> J/ip tttt) = (48.26 ± 0.95)%, (10) 



leads to 



r(^' ^ J/%Ij tttt) = 162.6 ± 9.4 keV. 



(11) 



aE can then be determined by comparing the value of 
(O with (HIl), and this leads to 



ue = 0.51. 



(12) 



The determination of is not so direct. So far the 
only experimental result can be used as an input da- 
tum to determine um /a e is the CLEO-c experiment on 
searching for the he state via ip' ^ hc + tt" [1J| . CLEO-c 
measured the product of the branching ratios [l^ 



B{il}' /ic7r°) X B{hc 



>rici) = (3.5 ±1.0 ±0.7) X 10""*. 
I 



(13) 



In an earlier paper [4| , we calculated the related branch- 
ing ratios by using QCD multipole expansion and the 
Gross- Treiman-Wilczek formula for obtaining the 



J 



transition rate T{^' — > hci:^), and by using perturbative 
QCD for obtaining the hadronic decay width V{hc — s- 
hadrons). Our obtained result is 4^ 5j 



B(V'' ^ hcTT°) X B{hc -> Vol) 



1.94 [ ^ 

aE 



X 10" 



(14) 



r 



Comparing these two results, we obtain 



as 



1.8 ±0.9. 



(15) 



Note that there is a 50% error in ([T5|) coming from the 
errors in (jl3p . 



We would like to mention that the small pion mass ef- 
fect [l^ is neglected in Eq. That the pion mass effect 
is negligibly small was argued in Rcf. 16]. In the present 
case, the uncertainty from the input value of am/ cue in 
Eq. is 50%. Compared with this large uncertainty, 
the small pion mass effect is evidently negilible in the 
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present calculation. 

With all these, we obtain the transition rate 

r(Xci ^ r/cTTTr) 11.0 (^^^ keV = 19.8 ± 9.9 keV.(16) 

In (fT6|) , the 50% uncertainty comes from the error in (fTSj) . 
For checking the model dependence of this prediction, we 
also calculated the transition rate using the well-known 
Cornell model with Coulomb plus linear potential (iTj . 
and the result is T{xci rjc'ni'n:) — 9A(aM I OiE) keV = 
17.0 ± 8.5 keV. So that the model dependence is about 
14% which is not significant considering the large un- 
certainty in The total width of Xci is r(xci) = 
0.89 ± 0.05 MeV [ll|. So that the branching ratio is 

B{xci JycTTTT) = (2.22 ± 1.24)%. (17) 

At e+e^ colliders, the state Xci can be produced via 
V'' ~^ IXci at the ip' peak. In the rest frame of -0', the 
momentum of Xci is 171 MeV which is only 5% of its 
mass M-^^^ = 3510.66 MeV. Therefore we can neglect the 
motion of Xci , and simply take the branching ratio ([T7]) 
to estimate the event numbers in the experiments. 

The detection of the process 

i>' IXcl JVcTTTT (18) 

can be performed in two ways, namely the inclusive and 
the exclusive detections [S]. In the inclusive detection, 
only the photon and the two pious are detected, while 
rjc is regarded as a missing energy. In the exclusive de- 
tection, all the photon, the two pions, and the decay 



products of rjc are detected. Reconstruction of 77c and 
Xci from the measured final state tagging particles can 
suppress the backgrounds. 

The inclusive detection requires measuring the mo- 
menta of the photon and the pions to certain precision. 
It is difficult to do this kind of analysis with the BES II 
data because the BES II photon detector is not efficient 
enough. At BES III and CLEO-c, this kind of detection 
is possible. 

For BES III, it is not difficult to accumulate 10* of tp' 
events. The branching ratio of ip' — > 7Xci is B^ip' — > 
jXci) = (8.7 ± 0.4)% [12. Taking account of a 15% 
detection efficiency, we obtain the number of events of 
(HHD at BES III 

BES III : 

N,nci{^' ^ 7Xci ^ JVcT^t:) = (2.90 ± 1.74) x 10^ (19) 

This is quite a large number. 

CLEO-c is now running at the ip' peak again, and will 
soon accumulate 3 x 10^ -0' events. The CLEO-c detec- 
tion efficiency is around 15% [l3|- For such a sample , 
the corresponding number of events of (|18p is 

CLEO - c : 

N,nci{i^' ^ 7Xci ^ IVcTttt) = (8.7 ± 5.2) x 10^. (20) 
So the transition p8)) can also be clearly studied at 
CLEO-c. 

For the exclusive detection, suitable decay modes of 
rjc should be taken for identifying the 77c in p^ . Some 
feasible decay modes with reasonable branching ratios 




Vc - 


pp ■■ 


B{Vc- 


^ pp) = (2.0 ±0.7)%, 


(21) 


Vc - 


-> K*K* : 


B{Vc- 


. K*K*) = (1.03 ±0.26)%, 


(22) 


Vc - 




BiVc- 


00) (0.27 ±0.09)%, 


(23) 


Vc ' 


^ K*{892)"K-n+ : 


B{Vc- 


K*{%92fK-T:+) = (2.0 ± 0.7)%, 


(24) 


Vc - 


K+R-TT+n- : 


B{Vc- 


K+R-TT+Ti-) = (1.5 ± 0.6)%, 


(25) 


Vc - 


2(7r+7r") : 


B{Vc- 


-> 2(7r+7r-) = (1.20 ±0.30)%, 


(26) 


Vc - 


■q'K'K — > ^JTTTT : 


B{Vc- 


T^TTTT 777r7r) = (1.9 ± 0.7)%, 


(27) 


Vc - 


KsK^n^ : 


B{Vc- 


-^KsK^TT^) = (1.9 ±0.5)%, 


(28) 


Vc ' 


-> KlK^tt^ : 


B{Vc - 


KlK^tt^) = (1.9 ± 0.5)%. 


(29) 



The modes ([281) and (HH) have been used by CLEO-c in 
the search for the he state [l3l |. 

In the exclusive detection, the requirement of the pre- 
cision of the photon momentum measurement in (|18p is 
not so strict. So it is possible to do this kind of analysis 
with the BES II data except for those modes with pho- 
tons in 77c decays. Therefore we also calculate the event 
numbers at BES II. The BES II data contains 1.4 x 10^ 



0^' events. Considering the ability of the BES II detector, 
we take a detection efficiency of 10% for BES II. The Pre- 
dicted event numbers for BES II, BES III, and CLEO-c 
are listed in TABLE II. We see that the exclusive detec- 
tion of the process (fT8|) can be well studied at BES III 
and CLEO-c. Considering the theoretical uncertainty, 
the events at BES II are marginal. 
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TABLE II: Predictions for the exclusive detection event numbers for the process (|18p in the CK potential model using the rjc 
decay modes shown in (I21|l — (|29[l at BES II, BES III, and CLEO-c. The accumulated numbers of the ip' events are taken to 
be 1.4 X 10'' for BES II, 10* for BES III, and 3 x 10^ for CLEO-c. The detection efficiency is taken to be 10% for BES II, and 
15% for BES III and CLEO-c. 



modes 


BES II 


BES III 


CLEO-c 


PP 


54 ±51 


579 ± 552 


174 ± 165 




28 ± 24 


298 ± 255 


90 ±76 


H 


7±7 


78 ±73 


23 ± 22 


K* (892)° K--K+ 


54 ±51 


579 ± 552 


174 ± 165 


K+K--K+TV- 


41 ±41 


435 ± 435 


130 ± 130 


2(7r+7r-) 


32 ± 28 


348 ± 296 


104 ± 89 


Tjirn 777r7r 


51 ±50 


550 ± 534 


165 ± 160 




51 ±44 


550 ± 476 


165 ± 143 




51 ±44 


550 ± 476 


165 ± 143 



We have seen that the theoretical uncertainty mainly 
conies from the experimental errors in (|13p . If the prod- 
uct of branching ratios B{ip' hcir^) x B{hc — > »7c7) can 
be measured to higher precision at CLEO-c or BES III, 
the predictions for the event numbers of the process p8| 
will be more definite. 

Finally, we would like to mention that the numbers 
listed in TABLE II are obtained by taking the PDG best 
fit experimental value r(i/'') = 337 ± 13 keV in Eq. ^ 
as input. Actually, the Particle Data Book also provides 
a world averaged value r(i/)') — 277 ± 22 keV in Meson 
Particle Listing [l^ . If we take this world averaged value 



as input instead of the PDG best fit value in ([TO]) , all the 
numbers in TABLE II will be smaller by a factor of 0.82. 
This is another uncertainty caused by the input datum 
of T{ip'). We expect further improved measurement of 
F (?/;') to reduce this kind of uncertainty. 
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